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ABSTRACT 

(N : 
o . 

O ■ We present an re-analysis of the longest timescale gravitational microlensing event discov- 

ered to date: MACHO-99-BLG-22/OGLE-1999-BUL-32, which was discovered by both the 
MACHO and OGLE microlensing alert systems. Our analysis of this microlensing parallax 
'"J . event includes a likelihood analysis of the lens position based upon a standard model of the 

0> I Galactic velocity distribution, and this impUes that the lens could be a black hole of ~ lOOM© 

I at a distance of a few hundred parsecs in the Galactic disk or a massive stellar remnant (black 

^ I hole or neutron star) in the Galactic bulge. Our new analysis includes data from the MACHO, 

\^ • GMAN, and MPS collaborations in addition to the OGLE data used in a previous analysis by 

O ■ Mao et al (2002). The crucial feature that distinguishes our analysis from that of Mao et al is 

. an accurate constraint on the direction of lens motion and an analysis of the impUcations of this 

O , direction. 
(N . 
O ■ 

: 

^! 1. Introduction 

o ■ 

^— > ■ The abundance of long timescale microlensing events towards the Galactic bulge has remained a puzzle 

• (Han & Gould 1996; Bennett et al. 2002) since shortly after the first systematic analyses of Galactic bulge 

^ ! microlensing events were reported (Udalski et al. 1994a; Alcock et al. 1997b). A long time scale event can be 

^ , caused by a massive lens, a slow relative transverse velocity between the lens and source, or by large source 

H ' star-lens and lens-observer separations. Thus, the excess of long timescale events could be caused by errors 

■ - - ' in our assumptions regarding the phase space distribution of the source stars or the lenses, or they could 

be caused by an unexpectedly large population of massive lenses. However, the phase space distribution 
of stellar mass objects in the Galactic disk and bulge is tightly constrained by observations, whereas little 
is known about the mass function of massive stellar remnants. Thus, it is possible that the excess of long 
timescale events is due to a population of black holes and neutron stars that is larger than expected. 

This possibility has recently received observational support from the MACHO (Bennett et al. 2002) and 
OGLE (Mao et al. 2002) collaborations which have identified individual candidate black hole microlensing 
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events based upon the detection of the microlensing parallax effect (Refsdal 1966; Gould 1992; Alcock et 

al. 1995). This refers to the detection of light curve features due to the Earth's accelerating motion around 
the Sun, and it is often detectable in long timescale microlensing events. To date, some 12 clear cases of this 
effect have been reported by the MACHO (Alcock et al. 1995; Bennett et al. 2002; Becker 2000), OGLE 
(Mao 1999; Soszynski et al. 2001; Smith, Mao, & Wozniak 2001; Mao et al. 2002), MOA (Bond et al. 
2001) and PLANET (An et al. 2002) collaborations. These microlensing parallax events are useful because 
they allow one to learn more about the nature of the lens than can be learned for most microlensing events. 
For most microlensing events, the only measurable parameter which constrains interesting properties of the 
lens is the Einstein diameter crossing time, t, which depends on the lens mass (M), distance (D^), and the 
magnitude of the transverse velocity {v± = \v±\). It is given by 



t = ^ = A. GMDe{Ds-De) 

where Dg refers to the distance to the source (typically ~ 8kpc for a bulge source), and Re is the radius of 
the Einstein Ring. In a microlensing parallax event, it is also possible to measure v, which is the lens star's 
transverse speed projected to the Solar position, given by 

V = VxD,/(Ds - D,) . (2) 

The measurement of both t and v give us two measurements for three unknowns (M, v±, and D^), which 
allows us to solve for the mass as a function of distance: 

16G DiDs 16G xDs ' 

where x = De/Ds. For the MACHO-99-BLG-22/OGLE-1999-BUL-32 microlensing event, (Mao et al. 
2002) have used Eq. 3 to argue that the lens mass is well above the measured mass of neutron stars 
I.4M0) for plausible values of x. 

Bennett et al. (2002) have also made use of the measured direction of the projected velocity, v, in 
a likelihood analysis employing the velocity and density distributions of a standard Galactic model. This 
allows a likelihood estimate for the lens distance, x, which can be converted to an estimate of the lens mass 
using Eq. 3. This analysis has revealed two other candidate black hole lenses, MACHO-96-BLG-5 and 
MACHO-98-BLG-6, as well as several other events with possible neutron star lenses. In this paper, we 
add additional data to the analysis of the MACHO-99-BLG-22/OGLE-1999-BUL-32 microlensing event, 
and show that this results in tighter constraints on the microlensing parallax parameters. We then apply the 
likelihood analysis developed in Alcock et al. (1995) and Bennett et al. (2002) to estimate the likely distance 
and mass of the lens. 
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2. The Data Set and Previous Analysis 

The MACHO-99-BLG-22/OGLE-1999-BUL-32 microlensing event was discovered and announced 
by the MACHO Alert system^ (Alcock et al. 1996), and it was independently discovered by the OGLE 
early warning system (Udalski et al. 1994b) about two months later. Following the MACHO alert, this 
event was added to the observing schedule for the GMAN (Becker 2000) and MPS (Rhie et al. 1999) mi- 
crolensing follow-up groups. MACHO and GMAN observations ended in 1999, while MPS and OGLE 
observations continued through 2000. Photometry from the images of this event is particularly challeng- 
ing because of severe crowding within ~ 1.5" of the source star. This makes it difficult for point spread 
function fitting, crowded field photometry programs such as DOPHOT (Schechter, Mateo, & Saha 1993), 
DAOPHOT/ALLFRAME (Stetson 1994), or SoDOPHOT (Bennett et al. 1993; Alcock et al. 1999) to obtain 
accurate photometry for this event. For the GMAN and MPS images, we were able to solve this difficulty 
by carefully selecting template images with sub-arc second seeing, but no such images are available for the 
MACHO data. There is a great deal of scatter in the MACHO-red data which appears to be inconsistent 
with the other data sets. A similar problem was seen in the MPS data, but it disappeared when the MPS 
data was reduced with a better seeing image as a template. This suggests that the photometric problem with 
the MACHO-red data is due to close neighbor stars that are not quite resolved in the MACHO-red template. 
The MACHO-blue data has much less scatter, and so we have decided to use only the MACHO-blue data in 
our analysis. 

The data set presented in this paper consists of 609 MACHO-Blue band measurements? 245 OGLE 
1-band measurements from their catalog of difference imaging measurements (Wozniak et al. 2001), 179 
MPS R-band observations from the Mt. Stromlo 1.9m telescope, and 21 GMAN R-band observations from 
the CTIO 0.9m telescope. The MPS data are reduced using a slightly modified version of the MACHO 
photometry code, SoDOPHOT. The photometric errors reported by SoDOPHOT are modified by adding a 
1.0% assumed systematic uncertainty in quadrature. The CTIO data were reduced with the ALLFRAME 
package (Stetson 1994), with the error estimates multiplied by a factor of 1.5 to account for systematic 
errors. 

A microlensing parallax fit for this event was previously published by Mao et al. (2002). However, Mao 
et al. (2002) present their fit results in a coordinate system that is inconvenient for Galactic bulge microlens- 
ing events. The difficulty is that all geometrical quantities are projected into the ecliptic plane, rather than 
a plane perpendicular to the line of site to the lens that is used for most discussions of microlensing events. 
Since the ecliptic plane intersects the Galactic bulge, this projection can be very nearly singular for Galactic 
bulge events. Instead, if we project the Earth's motion in the ecUptic plane to the plane perpendicular to 
the line of site (as in Alcock et al. (1995), for example), then there is no singularity. Event MACHO-99- 
BLG-22/OGLE-1999-BUL-32 is about 5 degrees from the ecliptic plane (A = 271.122°, /3 = -5.142°). 

'a catalog of MACHO alert events is available at http://darkstar.astro.washington.edu/ along with some 
information regarding eacli event. However, the reported timescales of these events are systematically underestimated because they 
come from fits that do not allow for blending. 

^The MACHO Project survey data are available from http://wwwmacho.mcmaster.ca/ and 
http : / / wwwmacho . anu . edu . au. 



-4- 



This means that the transformations of the error bars from the coordinate system of Mao et al. (2002) to our 

coordinate system are highly non-linear, and it means that the 1 a error bars reported in Mao et al. (2002) 
cannot be used to estimate the 2 a error bars. Because of this difficulty and the fact that Mao et al. (2002) 
use a data set which differs from the publicly available one, we will not attempt a detailed comparison of 
our fit parameters with their results. 



3. Microlensing Parallax Analysis 

The magnification for a normal microlensing event with no detectable microlensing parallax is given 

by 

«2 + 2 



A(t) = ■ u{t) = ^ul + [2{t - t,)/t)Y , (4) 

where to is the time of closest approach between the angular positions of the source and lens, and uq = b/ Re 
where b is the distance of the closest approach of the lens to the observer-source line. Eq. (4) can be 
generalized to the microlensing parallax case (Alcock et al. 1995) by assuming the perspective of an observer 
located at the Sun. We can then replace the expression for u{t) with 

u^(t) = ul + uj^{t -tof + SiTl'[Q.{t -tc)] 

+2a sin[r2(t — tc)] [oj{t — to) sin 6 + uq cos 9] 

+a^sm^ 13 cos^[Q,{t - tc)] + 2a sin /3 cos[r2(t - tc)] [uj{t - to) cos 6 - ■uosin6'] (5) 

where A and (3 are the ecliptic longitude and latitude, respectively, 6 is the angle between v± and the North 
ecliptic axis, lo = 2/t, and tc is the time when the Earth is closest to the Sun-source line. The parameters a 
and $7 are given by 

a = ^^^^^ (1 - ecos[J^o(t - t^)]) , (6) 

and 

n{t - tc) = no{t - tc) + 2esin[no(t - tp)] , (7) 

where tp is the time of perihelion, Uq = 27r yr~^, e = 0.017 is the Earth's orbital eccentricity. 

The microlensing parallax fits were performed with the MINUIT routine from the CERN Library, and 
the results are summarized in Table 1, and are shown in Fig. 1 The first four fit parameters listed in Table 1 

are fit blend fractions for each of the four data sets. Because the OGLE data was reduced with a difference 
imaging method, there is no photometric reference flux value to use to define /ogle = 1- Therefore, we 
select the fit flux for the best parallax fit to define /ogle = 1- 

Although the comparison with the parallax fit of (Mao et al. 2002) is complicated by their peculiar 
choice of coordinates, we can compare to their t and v values since these are not affected by their coordinate 
system. They find t = 1280 days and v = 79 ± 16km/s, while we find t = 1120 ± 90 days and 
V = 75^g km/s. We find that the parallax fit gives an overall improvement of 541.7 which compares 
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to the improvement of 298.1 found by (Mao et al. 2002). Thus, our results seem to be somewhat more 
accurate, but consistent with the Mao et al fit. 

The microlensing parallax fit parameters can be used to determine the lens mass as a function of dis- 
tance following Eq. 3 as shown in Fig. 2. 

3.1. Likelihood Distance and Mass Estimates 

We can make use of our knowledge of the phase space distribution of lens objects in the Galactic 
disk and bulge to further constrain the distance and mass of the lens because the probability to obtain the 
measured value of v is quite sensitive to the distance to the lens. Following Alcock et al. (1995) and Bennett 
et al. (2002), we can define a likelihood function 

L{x; v) oc ^/x{1-x)pl{x) v{1 - xf J fs{vs) /l((1 - x){-Vq + v) + xvs) dvs, (8) 

where is the density of lenses at distance x = D^/Ds, and the integral is over combinations of source 
and lens velocities giving the observed v. V5 and = (1 — x) (v0 + v) + XV5 are the 2-D source and lens 
velocity distribution functions. We assume the following Galactic parameters: a disk velocity dispersion 

of 30km/s in each direction, a flat disk rotation curve of 200km/s, and a bulge velocity dispersion of 
80km/s with no bulge rotation. The density profiles are a standard double-exponential disk and a Han & 
Gould (1996) barred bulge, and the source star is assumed to reside the the Galactic bulge at a distance of 
8.5 kpc which is consistent with the measured radial velocity, Vr = -1-38 it 2 km/sec of the source star (Cook 
et al. 2002). 

The resulting UkeUhood functions for is shown as the long-dashed curves in Fig. 2. This differs 
from the UkeUhood functions for the six events presented by Bennett et al. (2002) because it has two peaks. 
This is because the measured direction of v is nearly opposite of the direction of Galactic disk rotation. For 

a typical Galactic disk lens, we expect v to be roughly parallel to the disk rotation direction as is observed 
for the six MACHO events analyzed in Bennett et al. (2002). Since the bulge is assumed to have little or 
no rotation, the average Une-of-site to a bulge lens is stationary at the position of the source. (A rotating 
bar-shaped bulge could imply an average bulk motion at along some Unes-of-sight.) The flat rotation curve 
of the Galactic disk will cause stars orbiting interior to the Solar circle to have a higher angular velocity 
than the Sun, and so the preferred projected velocity, v will be in the direction of Galactic disk rotation. 
However, the probability for v to be in the anti-preferred direction is not small (Bennett et al. 2002), so the 
fit parameters for MACHO-99-BLG-22/OGLE-1999-BUL-32 should not be considered unlikely 

The double peaked nature of the likelihood function shown in Fig. 2 can easily be understood if we 
consider the ways in which a v parallel to the Galactic disk rotation can be avoided. One possibility is for 
a disk lens to be very close to us. In this case, the average Galactic angular velocity at the lens position is 
only slightly larger than that of the Sun, and so the velocity dispersion of the disk wiU give a substantial 
population of potential lens objects with retrograde v values. The other possibiUty is that the lens resides 
in the Galactic bulge where there is Uttle or no systematic rotation. In this case, with both the lens and 
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source in the bulge, the Galactic orbital motion of the Sun gives a preference for a retrograde v value, but 
the preference is fairly small because the lens must be much closer to the source than to the Sun. The two 
likelihood function peaks can be then understood as being due to disk and bulge lenses. This is clear from 
the decomposition of the likelihood function shown in Fig. 2. For a disk lens, the likelihood function peaks 
at a mass of 130Mq, which is about an order of magnitude above the most massive of the known stellar 
mass black holes in x-ray binary systems. If the lens is in the bulge, however, the most likely mass is only 
about 4Mq. 

One common way to interpret likelihood functions is the Bayesian method, in which the lens mass 
(or distance) probability distribution is given by the likelihood function times a prior distribution, which 
represents our prior knowledge of the probability distribution. In our case, the likelihood function represents 
all of our knowledge about the lens mass and location, so we select a uniform prior With a uniform prior, 
the likelihood function becomes the probability distribution, and we are able to calculate the lens mass 
confidence levels listed in Table 2. Also included in this table are the lens mass confidence levels if the lens 
is known to be a disk or a bulge star. If the lens resides in the disk. Table 2 indicates that the lens mass 
should be in the range I6M0 < M < 172M0 at 1 a and 6.3 < M < 9OOM0 at 2 a. This wide range is 
due to the fact that the distance to a disk lens is not well constrained. If the lens resides in the bulge, then 
the mass range is much narrower: 2.2M0 < M < 5.5M0 at 1 cj and 1.34M0 < M < 8.3M0 at 2 a. This 
suggests a low mass black hole, although a neutron star or even a main sequence lens cannot be ruled out. 

3.2. Constraint on a Possible Main Sequence Lens 

If the lens star is a main sequence star, then there is an additional constraint on the distance and mass 
of the lens due to the observed upper limit on the brightness of the lens star (Bennett et al. 2002). We have 
assigned a conservative upper limit on the V-band brightness of F > 20.1 it 0.5 based upon the best parallax 
fit. We have applied this constraint to the likelihood function by multiplying the likelihood function by the 
Gaussian probability that the lens brightness exceeds upper limit on the lens star brightness. The result is 
shown as the alternating short and long-dashed curve in Fig. 2, which indicates that main sequence lens stars 
more massive than about I.5M0 is ruled out. A much tighter constraint on the brightness of the lens can 
be obtained with HST images as has been done with the MACHO-96-BLG-5 black hole candidate (Bennett 
et al. 2002). 

4. Discussion and Conclusions 

We have presented a microlensing parallax analysis of the MACHO-99-BLG-22/OGLE-1999-BUL-32 
microlensing event using the publicly available OGLE and MACHO data along with MPS and GMAN data. 
The projected velocity indicated by the microlensing parallax fit has a direction that is nearly opposite of 
the direction of disk rotation. This suggests that the lens is either very close to us in the Galactic disk or in 
the Galactic bulge. The microlensing parallax mass-distance relationship, 3, indicates that the lens must be 
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very massive if it is nearby, with a mass of order ~ lOOM©. If the lens is in the Galactic bulge, however, its 
mass is much smaller, ~ 4M0. In both cases, the lens is likely to be a black hole. These conclusions have 
been quantified with a likelihood analysis similar to those presented by Alcock et al. (1995) and Bennett 
et al. (2002). We should note that this analysis does depend on the assumption that stellar mass black holes 
are not bom with a large velocity "kick'.' This assumption is supported by the available evidence (Nelemans, 
Tauris, & van den Heuvel 1999), but if this assumption is wrong, then the MACHO-99-BLG-22/OGLE- 
1999-BUL-32 could be a ~ lOM© black hole half way to the Galactic with a rotation speed of only about 
60 km/s. 

We have not carried a Bayesian likelihood analysis with an assumed mass distribution prior as has been 
advocated by Agol et al. (2002) and also done by Bennett et al. (2002). Such an analysis is of little use for 
this event because the lens mass is very likely to be > 1.5Mq. Since a main sequence lens of this mass is 
excluded, we would have to provide a prior distribution of stellar remnants, but there is virtually nothing 
known about mass function of stellar remnants above l.SM©. 

Our analysis agrees with the previous analysis of Mao et al. (2002) and indicates that the MACHO- 
99-BLG-22/OGLE-1999-BUL-32 lens is very likely to be a black hole, but we show that a nearby, massive 
black hole of ~ lOOM© and a low mass black hole of ~ 4M0 in the Galactic bulge are two distinct 
possibilities that can explain the observed microlensing parallax parameters. When this event is combined 
with the two other candidate black hole microlenses presented by the MACHO Collaboration (Bennett et al. 
2002), it appears that black holes may comprise a significant fraction of the mass of the Galactic disk and 
bulge, perhaps as large as ~ 10%. Future observations with large ground based interferometers (Delplancke, 
Gorski, & Richichi 2001) should be able to accurately determine the mass of future black hole microlenses 
and determine the black hole mass fraction. 

This work was supported, in part, through the NASA Origins Program Grant NAG5-4573. 
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Fig. 1.— The light curve of the MACHO-99-BLG-22/OGLE-1999-BUL-32 microlensing event with data 
from MACHO, OGLE, MPS and GMAN. The blue, green, red, and magenta data points indicate data from 
the MACHO-Blue band, the Mt. Stromlo 1.9m R-band,the CTIO 0.9m R-band, and the OGLE I-band, 
respectively. The solid curve shows the best fit microlensing parallax event while the cyan curve shows the 
best fit standard microlensing light curve. In the upper panel, the data have been averaged into 4-day bins, 
while in the lower panel the data have been averaged into 10-day bins. 
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Fig. 2.— The mass vs. distance relation (solid curves) for the MACHO-99-BLG-22/OGLE-1999-BUL-32 
microlensing event is shown along with the likelihood functions (long dashed curves) computed assuming a 
standard model for the Galactic phase space distribution. The (green) short dashed curve is the contribution 
from disk lenses, and the (red) dotted curve is the contribution from bulge lenses. The alternating long 
and short dashed curve is the likelihood function which includes the constraint on the brightness of a main 
sequence lens star. A main sequence lens is quite unlikely for this event. 



- 11 - 



Table 1 . Microlensing Fit Parameters 



Fit 


/ogle 


/mb 


/CTIO 


/mps 


to (MJD) 


wo 


t (days) 


V (km/sec) 


^ (dof) 




Par. 


1.00(11) 


0.26(3) 


0.44(5) 


0.46(5) 


2745.0(1.5) 


0.043(6) 


1120(90) 


75(9) 


1.85(20) 0.864 


899.8 


Std. 


0.02(5) 


0.005(3) 


0.01(1) 


0.01(1) 


2742.4(1) 


0.0007(1) 


4.0(2) X 10^ 


= 


= 1.382 


1441.5 



Table 2. Microlensing Parallax Likelihood Mass Estimates 



lens location 


P = 2.5% 


P = 5% 


Confidence Levels P(M/M0 < N) 
P = 16% P = 50% P = 84% 


P = 95% 


P = 97.5% 


disk or bulge 


560 


310 


106 


17.1 


3.2 


2.1 


1.70 


disk 


900 


480 


172 


49 


16.2 


8.5 


6.3 


bulge 


8.3 


7.3 


5.5 


3.6 


2.2 


1.60 


1.34 



